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Gauge conditions for moving black holes without excision 
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Recent demonstrations of unexcised, puncture black holes traversing freely across computational 
grids represent a significant advance in numerical relativity. Stable and accurate simulations of 
multiple orbits, and their radiated waves, result. This capability is critically undergirded by a careful 
choice of gauge. Here we present analytic considerations which suggest certain gauge choices, and 
numerically demonstrate their efficacy in evolving a single moving puncture. 

PACS numbers: . 


I. INTRODUCTION 


Gravitational waveforms from binary black hole merg- 
ers can only be obtained through 3D numerical relativity 
simulations of the full Einstein equations. Such simula- 
tions have proven to be very challenging. Improvements 
have often come from new formulations or a new choice 
of gauge. For example, a corotating gauge led to the first 
evolution of a binary through one full orbit [1] as well as 
through a plunge, merger, and ringdown [2]. 

Recently, the authors of this paper [3] and another re- 
search group [4] independently developed the capability 
to numerically evolve freely moving black holes, Le. black 
holes that are in no way fixed to a any position in the 
computational grid, without using excision techniques for 
the black hole interiors. These efforts led to the simula- 
tion of record breaking numbers of orbits, among other 
accomplishments. Both of these research programs em- 
ployed gauge conditions that differed not only from those 
prescribed in prior literature, but also from each other. 
A careful choice of gauge is important, particularly for 
moving black holes, as it can determine whether a sim- 
ulation is stable or becomes irrevocably corrupted by a 
runaway build up of numerical error. There are various 
ways in which a choice of gauge can affect the stability 
and accuracy of a numerical simulation; for example, a 
poor choice of gauge might cause large gradients in the 
fields which in turn engender large finite differencing er- 
ror. 

In this paper we investigate a group of gauge choices 
analytically, and numerically demonstrate their effects on 
moving black hole simulations. These gauge choices are 
framed in the context of the conventional BSSN formu- 
lation of Einstein’s equations. After a brief exploration 
of the shortcomings of a commonly used gauge choice in 
II, we calculate in section III A the characteristic speeds 
for various gauges and then show the numerical behavior 
of some of them in section III B. In section IV we discuss 
the properties of the most promising gauge choices and 
suggest further improvements. Conclusions are presented 
in section V. 


II. ORIGINAL 1+LOG SLICING AND 
GAMMA-DRIVER SHIFT CONDITIONS 


We consider first a gauge introduced in [? ' 
ing non-moving punctures without excision: 

| for evolv- 

d t a = -2 aK 

(i) 

dtP 1 = \a% n B i 

(2) 

d t B i = d t r - r j B i 

(3) 


where a is the lapse, fi l is the shift, K and f* are the 
usual BSSN variables, is the initial conformal factor 
P ])> IS a . n auxiliary variable to make the shift 
equation first order in time, and r] is a constant damp- 
ing factor. This gauge has proven to be very effective 
at evolving a single non-moving puncture for arbitrarily 
long durations, and in some circumstances it can also be 
applied to binary punctures. Here the 14-log condition 
on the lapse is a standard singularity- avoiding slicing 
condition. Of particular interest is the “Gamma-driver” 
condition for the shift, which is designed to ultimately 
“freeze Gamma”, i.e. drive d t f to zero as the physics of 
the spacetime also evolve towards quiescence. 

For n > 0, n = 0 at the puncture, and thus this spe- 
cific gauge will ensure that fi % will not evolve there. As 
the motion of the puncture is determined by the shift, the 
puncture will remain motionless in this gauge since the 
shift is initially zero. For black holes that are physically 
dynamical the coordinates may become highly distorted. 
For example, in the case of binary black holes, the separa- 
tion between the black holes as determined by the proper 
distance from horizon to horizon in some spatial slice 
must decrease with time. This implies that components 
of the 3-metric will approach' zero between the black 
holes. Then the evolved part of the BSSN conformal 
exponent <fi must either approach — oo, or components of 
the conformal 3-metric j \j must approach zero. In the 
latter case the inverse 3-metric must diverge in order to 
maintain a unit determinant. Thus, in the fixed puncture 
approach, arbitrarily large, ever-increasing fields and gra- 
dients around the black holes seem likely to develop and 
cause numerical difficulties. Indeed this can be demon- 
strated even in the case of a single moving black hole. 
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FIG. 1: Results of the non-moving puncture gauge dtot = 

~2aK , dtff- = §a^“ n £\ = F - at time * = 17M. 

The conformal metric becomes singular by t — 18 M. 

We will use as a test case a single puncture black hole, 
as given by Bowen- York initial data, with unit puncture 
mass and momentum such that its physical speed should 
be half the speed of light. The puncture position will 
start at coordinates (—3 M, 0, 0) and the momentum will 
be in the ^-direction. and B l will be initialized at 
zero, while the lapse will be precollapsed with the profile 
2 as suggested in [? ] (and recommended in f? ] for 
moving punctures). The damping parameter rj will be 2, 
unless otherwise stated. We will assess the performance 
of each gauge by the evolving behavior of the quantities 
a, /?*, 0, and T x ; in particular the peak in <p will roughly 
indicate the position of the puncture and any extreme 
gradients in F x ~ —dj^ x will tend to indicate problems 
with the gauge. 

In this regard it should be noted that F is completely 
controllable by the gauge condition, and so any unde- 
sirable features observed in F are in principle avoidable 
via a better choice of gauge. In particular, the specific 
“Gamma-freezing” condition f 1 — 0 may be desirable, as 
this elliptic generalization of isotropic coordinates pro- 
posed by Dirac is expected to be non-pathological. So, 
among hyperbolic Gamma-driving conditions, those that 
result in smaller values of F and its derivatives might be 
preferred. 

In Fig. 1, evolving with the original gauge given in 
Eqs. (3), we see that T x grows very large and (p devel- 
ops sharp features. Shortly thereafter the inverse confor- 
mal metric diverges, which is our criterion for stopping 
the run. This failure motivates investigation of moving 
punctures. 

If we wish the punctures to move we must allow 
to evolve freely away from zero and thus we must take 
the conformal factor out of the shift equation (setting 
n = 0 in Eq. (2). The results, shown in Fig. 2, are not 
much improved. The puncture is' now free to move but 
non-propagating features at x = —3 M destabilize the 


FIG. 2: Results of the gauge dta = -2aK, dtp * = f aB\ 

dtB 1 ±= F — rjB 1 at time t = 17 M. Sharp features around 
x = —3 M fail to propagate. The conformal metric becomes 
singular by t == 18 M. 


simulation. Clearly there is a problem with zero-speed 
modes. The extended, collapsed region of the lapse points 
to a particular difficulty in this regard, and suggests an 
improvement. 

The factor of a in Eq. (2) should also be removed, for 
even if it were not precollapsed it would still be driven 
nearly to zero around the puncture by the l+.log slic- 
ing condition, and therefore may retard the evolution of 
P z . When this factor is removed, the results are not as 
catastrophic as before, but the simulation is plagued with 
noise. Fig. 3 shows, in particular, a tendency for the min- 
imum in the lapse to lag behind the puncture position. 
The above results suggest that a careful study of the 
propagation speeds in this system might be helpful. 


III. CHARACTERISTIC SPEEDS 
A. Eigenvalue analysis 

A potential danger of zero-speed modes is that errors 
may couple to them, compound in place, and become 
inordinantly large. Zero-speed modes seem particularly 
hazardous in the case of moving black holes as they may 
adversely affect the dynamics by putting a drag on the 
puncture motion. In some cases zero-speed modes can 
also be a source of instability ([? ]). So there is good 
reason to be aware of such modes. 

For the purpose of investigating characteristic speeds, 
a simple yet revealing method of analysis is to linearize 
the equations, assume plane wave solutions, and solve for 
the eigenmodes of the resulting algebraic system. With 
the gauge conditions described above, the BSSN equa- 
tions can be linearized about flat space as follows. 





where 




FIG. 3: Results of the gauge dta = — 2 aK, dtf? = §£\ 

dtB* ■■= P — 77 IT at times £ = 4.5M (top panel) and t — 30 M 
(bottom panel). The top panel clearly shows the lapse lagging 
behind the puncture, while the bottom panel exhibits noise 
propagating from the puncture. 


a -- 

B { = 
b l - 

* = 
k = 

hij = 



—2aoK (4) 

r (5) 

( 6 ) 

(7) 

-didia + fiodkK (8) 


—2aoAij 4- diV 4- djb 1 — gbijdkb k + Pq dkhij(9) 
-didjd - la 0 dkdkhij + ~a 0 dif j 4- ^a 0 dj ,f* 


-2a 0 di<9j<£] TF 4- ffidkAij (10) 

f* = -\oi Q diK + dAV + + ffidkF ( 11 ) 


where a = a — ao, b l ~ /?* — /?g, hij = 7 ^ — 5^ and 
the damping term on B l has been dropped as it should 
not affect the real part of the characteristic speeds. The 
lapse and shift are assumed to have constant, zeroth order 
terms ao and /?q, respectively which conveniently allows 
us to capture essential effects of the lapse factors and the 
advection terms without making the eigenvalue problem 
intractable. 

For spatial variation only in one dimension, with no 
initial transverse components, and assuming plane-wave 
solutions, the above system of equations can be written 
in the form: 


( 12 ) 


l«> 


a V 
B 

b 


K 

h 


* 


gi(kx-cjt) 


(13) 


with a, B, 5, </>, JT, h 7 A, and f the amplitudes of a, B x , 
b x , (j>, h XX ) A xx , and t x respectively, and M == 

0 0 0 0 -|a 0 0 0 0 \ 

0 ^ 0 § 0 0 -Pq 

0 0 0 0 0 0 0 

0 0 -fa -efeao 000 

-ik 0 0 0 -Po 0 0 0 

0 0 -f 0 0 —Pq -%a 0 0 

■ _2|fe o 0 -f^ 0 0 -fa 0 -p 0 -|a 0 

0 0 ^ 0 fa 0 0 0 -po J 

(14) 

The eigenval ues of M are 0,— /3 q, — Po + a p,— /3p — 
a o>” 2 ^° 2 V $0 + 8ao,“|/?o •“ I V /?o + 8ao ? — §/?o + 

+ 4 a o +1 > and ~§A> ~ \\JPl + Aa Q +l - If («} is the 
eigenvector associated with eigenvalue u, and the covec- 
tor (?;| is defined such that J2 V = K then the eigen 
decomposition of M allows the system of evolution equa- 
tions to be written as a series of advection terms, each 
associated with a characteristic velocity equal to one of 
the eigenvalues: 


®t\u) = y^v9 x |t;)(u|w) (15) 

The inner product with (t>| gives: 

d t (v\u) = vdx(v\u) (16) 

Now we are interested in slow speed modes. Note 
that although ft 1 is typically initialized at zero, with the 
Gamma-driver condition it evolves rapidly and signifi- 
cantly enough from zero that we do not classify fto as a 
slow mode (since we’re interested in moving punctures, 
fi l is generally non- vanishing). In addition to the obvious 
zero speed mode there are two modes that approach zero 
in the limit where ao approaches zero. If |ao| << |/?o| 

then + ■+ 8a o and - §/? 0 + \ \/ Pi + 4ag +1 be- 

come approximately 4ao//?o and 2 aQ + 1 //?o, respectively. 
In the full nonlinear system, the 1+log slicing condition 
will collapse a nearly to zero in a finite region around the 
puncture (where /?* will be nonzero for a moving punc- 
ture). Modes which propagate at a speed approximately 
proportional to a positive power of a for small a, which 
we will refer to as a-speed modes, are thus of potential 
concern. In particular, late in the evolution of black hole 


d t \u) = ikM\u) 
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binaries there will emerge a significant, stationary region 
in which the lapse has collapsed nearly to zero, in which 
a-speed modes are effectively zero-speed m odes. 

The eigencovector for v = — |A) + + 8ao in- 

volves only the fields a and K 1 is relatively independent 
from the other fields, and easy to fix, so we will address 
it first. This a-speed mode can be understood also in the 
context of the original nonlinear equations by noting that 
a couples with K to give a wave equation such that the 
speed of fluctuations of the lapse goes to zero when the 
lapse itself goes to zero. It is particularly egregious when 
the lapse couples to an a-speed mode, as once it collapses 
nearly to zero it will tend to get “stuck” there. But this 
can be remedied in a natural way simply by adding an 
advection term as follows, 

d = -2 aK + I3 j dja, (17) 

which will modify the principal part of its wave equation 
such that when a goes to zero its sp eed will go t o /?. 

The eigencovector for -|A) +'§^£0 4-4a£ +1 is more 
complicated, involving not only a but several other fields 
as well. It turns out that advecting the lapse will not re- 
move this mode. But noting that this eigenfield involves 
b x , and that the equation for b x is conspicuously absent 
an advection term, an obvious strategem is to advect the 
original as follows: 


ft = Z -B l + ftdjft (18) 

Finally, Eq. (16) for the eigenvalue v = 0 yields 

d t (B x -T x ) ~ 0 (19) 

This equation, a trivial consequence of the original 
Gamma-driver condition, immediately suggests a simple 
modification. The most natural way to remove this zero- 
speed mode is to advect it at /? speed by modifying the 
B % equation thus: 

d t B i - d t r + P j d j (B i - f *) - rjB * (20) 

It turns out that the +(3 j djB t term alone will not suf- 
fice, as for example it might introduce an exponentially 
growing mode. Meanwhile the —/3 j djf z term alone adds 
another a-speed mode to the system. 

The foregoing analysis suggests the addition of four 
advection terms to the gauge equations, as indicated in 
Eqs. (17,??,??). The lapse equation is relatively inde- 
pendent of the others: if it is advected then a particu- 
lar a-speed mode is removed, and if it is not advected 
then this particular a-speed mode remains. But the 
three advection terms suggested for the shift equations 
are strongly interdependent, so their combined effects are 
not so obvious and it is instructive to consider which un- 
desirable modes result from which combinations of these 
three terms. The possible combinations and their result- 
ing “bad” speeds are summarized in Table 1. Here we 


see that only when all the above advection terms are 
included will the system be free of slow-speed and expo- 
nentially growing modes. In this case the benign eigen- 


speeds are: -Ay, An -/?o-ao ? “Ar+ao, 





— /?o — v/2oo, — A) + V2ao* 


case # 

+ftd j B i 

-ftdjt 1 

+P j d : P i 

“bad” speeds 

1 

N 

N 

N 

0, a p+1 

2 

N 

N 

Y 

0 

3 

N 

Y 

N 

a p+1 

4 

N ■ 

Y 

Y 

a p+1 

5 

Y 

N 

N 

0 

6 

Y 

N 

Y 

—i 

7 

Y 

Y 

N 

a p+1 

8 

Y 

Y 

Y 

none 


TABLE I: Effect of various advection terms in the shift equa- 
tions on the presence of undesirable eigenspeeds. ’Y 5 or ’N 5 
indicates whether each advection term is added or not, re- 
spectively, to Eq. (18) or Eq. (20). 0, a p+1 , or —i indi- 

cates whether the resulting linearized equations have a zero- 
speed mode, an a-speed mode, or an exponentially growing 
mode, respectively. The lapse is assumed to be advected as in 
Eq. (??); otherwise an additional a-speed mode would appear 
in every case. 


B. Numerical tests 


Evidence of the slow-speed modes found in the lin- 
earized analysis can often (but not always) be found, in 
simulations of the full nonlinear system. Note that zero- 
speed modes only represent a potential danger. In some 
cases errors may not couple to these modes even though 
they exist in the equations. 

Here we provide a few examples of the cases given in 
Table I. Fig. 4 represents case #1 with p — 1, which is 
predicted to have a zero-speed mode and indeed mani- 
fests a non-propagating feature at x = —3M. On the 
other hand, case #1 with p — 0, also predicted to have 
a zero-speed mode as well as an a-speed mode, belies 
no indication of either in Fig. 5, instead appearing quite 
smooth. In this case, evidently, the zero-speed modes do 
not couple significantly with the other fields. In Fig. 6 
the spike in f®, which becomes steeper in time as the 
“tail” in f* grows, may be related to the a-speed mode 
predicted for case #3. For a longer evolution this grow- 
ing gradient in f* can be expected^ to adversely affect 
the constraints. In Fig. 7 the tail in T x and in particular 
the non-propagating bend in the tail around x == —2 M 
seems to be evidence of the zero-speed mode predicted 
for case #5. In Fig. 8, noise in f® grows exponentially, 
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FIG. 4: Results for the gauge d t a = —2 aK 4- ftdjOt, dtfi x = 
f aB\ dtB 1 — dtF x - r]B x (case #1 in Table I with p=l) at 
time t = 14 M. Non-propagating features are evident around 
x = ~3M. The conformal metric becomes singular by t — 
15 M 



FIG. 5: Results for the gauge d t a = — 2 aK + fPdja, dtj? — 
| B % , d t B x = dtf x — 7]B x (case #1 in Table I with p=0) at 
time t = 30 M. The evolution appears smooth. 



x 


FIG. 6: Results for the gauge d t a — —2 aK -f ftdja, d t /3 x 
| oB% d t B i = d t r - FdjT 1 - rjB 1 (case #3 in Table I with 
p=0) at time t = 301^. F x continues to grow larger and its 
gradient steeper. 



x 

FIG. 7: Results for the gauge d t a ■== —2a K + dj a, dt/3 x = 
| aB x , dtB 1 = d t f 1 4- P j djB x — rjB x (case #5 in Table I with 
p~0) at time t — 30 M. The feature around x — —2 M does 
not propagate. 


as expected for case #6. And finally, case #8, depicted 
in Fig. 9 and Fig. 10, is demonstrated to be free of slow- 
speed modes .whether p = 0 or p ~ 1. (Although the for- 
mer appears much smoother, perhaps because the eigen- 
speeds -fio -j- y/oLQ allow faster propagation of error away 
from the puncture than — /?© 4- .a©-) 


IV. PROPERTIES OF THE “CLEANEST” 
GAUGES 

Of all the combinations suggested in Table I, two dis- 
tinguish themselves in numerical tests as conducive to 
particularly smooth propagation of the moving black 


hole. These gauges are case #8 with p = 0, which we 
will refer to as the “shifting-shift case” : 


dtF 

-0 j d j 0 i ~ 

\ Bi 


(21) 

dtB * • 

- 0 j d j B i = 

a,r- 

fidjF - rjB 1 

(22) 

and case #1 with p = 0, 

which 

we will refer to 

as the 

“non-shifting- 

shift” case: 





dtF = 

i b ‘ 


(23) 


d t B i = 

d,t‘- 

T)B i 

(24) 
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FIG. 8: Results for the gauge dta = —2 aK + p j dja , dtp 1 = 

+ p j d j B i - (case #6 in 
Table I with p=0) at time t = 10M The noise in r* grows 
exponentially. The conformal metric becomes singular by t == 
11M. 



x 


FIG. 9: Results for the gauge d t a — -2 aK -f P j dja, dtp 1 ~ 
faB { -F P j djP\ d t B 1 = dtf* + P j dj{B i - P) - 77, B 4 (case #8 
in Table I with p—1) at time t = 30M. Aside from the sharp 
features in P, which do not grow significantly in time, the 
propagation is non-pathological. 


The efficacy of Eqs. 23 in evolving moving punctures 
was originally demonstrated in [? ]. Apparently the zero- 
speed and a-speed modes are not problematical in this 
case. Regarding the second shift condition, Eqs. 21, we 
have recently used it successfully in the stable evolution 
of two equal mass black holes through 5.5 orbits. This 
latter gauge has also been recommended recently for the 
strong hyperbolicity it brings to the BSSN equations [? ], 
suggesting it may be a more robustly stable choice than 
Eqs. 23. Nevertheless for the single black hole simulations 
represented by Figs. 5 and 10, the results from these two 
gauge choices are nearly indistinguishable. 



FIG. 10: Results for the gauge d t a = -2aK + P j dja? dtp 1 = 
f B i + ftdjp 1 , dtB 1 = dtV -I- p j dj(B i - P) - r\B l (case #8 
in Table I with p=0) at time t == 30 M. The evolution is very 
smooth. 

For these two gauge choices, we now consider varying 
the damping parameter, which has up to now been set 
to T) = 2. By moving toward 77 — F 0 Eqs. (21) appear a 
little closer to realizing a “Gamma-freezing” condition. 
Indeed, with 77 = 0, we find that the the black holes 
come closer to realizing the physically expected velocity. 
Both the “shifting-shift” and “non-shifting-shift” gauge 
options with 77 = 0 have proven to allow stable evolutions 
of a single black hole. In Fig. 11 we find that in both cases 
setting 77 ~ 0 results in significantly smaller values for 
f*, meaning a closer approximation to the Dirac gauge. 
Further, we comment that we have observed numerically 
that for 77 > 0, f 4 appears to show very slow linear growth 
in time, whereas for 77 == 0, f z appears to be bounded. 
Comparing the “shifting-shift” and “non-shifting-shift” 
options with 77 = 0 we now find more noticeable differ- 
ences in P, with our recommended “shifting-shift” op- 
tion giving a smoother result with f 1 holding closer to 
zero near the puncture. 

These two gauges share an additional feature in com- 
mon, which is that either Eqs, (23) or Eqs. (21) can be 
integrated to give the relation 

B i = t i - ^ (25) 

since B % ~ F 1 = P % — 0 initially. This 'fact, also evi- 
dent numerically (Pig. 12), suggests substituting for B % 
in the evolution equation for /?*, to obtain (in the case of 
Eqs. (21)) 

= ^r + fidjP* -vP (26) 

Fig. 13 demonstrates that the resulting shift condition, 
Eq. (26), can be evolved numerically to yield the same 
stable, smooth simulation as with the analytically equiv- 
alent condition, Eqs. (21), used previously. This should 






7 



x 

FIG. 11: T x for the gauge d t a — —2 aK + ftdja, d t ft ■■== § B l : 
d t B z = d t f l (case #1 with r\ — 0) and for the gauge dta = 
-2aK+ftdja, ftjS 4 = f B'+ftdjft, d t B i = 
f ' l ) (case #8 with rj = 0) at time £ = 10M. In the “shifting- 
shift” case (case #8) T x is smaller and smoother around the 
puncture than in the “non-shifting-shift” case (case #1). 



x 

FIG. 12: Results for B x and f* — with the gauge dta = 
-2 aK+ftdja, d t ft = f B i +ftd j ft, d t B i = d t r + ft d^B 1 - 
f l ) — rjB * at time t = 15 M. The difference between the two 
curves is within 10~ n . 

hold true without the ft dj ft term. Dispensing with the 
B 2 evolution equation results in a slightly more efficient 
numerical implementation, and also guarantees removal 
of the zero-speed mode associated with 23\ 

Finally it may be noted that when rj — 0 and f * 0, 

Eq. (26) (and Eqs. (21)) admits a “shock” solution of 
the form (x ~ x 0 ) / (to — t) . This represents the nonlinear 
tendency of this equation to advect large magnitudes of 
the shift faster than smaller magnitudes, which for neg- 
ative values and a positive slope can lead to a vertical 
slope. However, as d t t 2 depends on derivatives of ft, T 2 
is unlikely to vanish when those derivatives become large. 



x 

FIG. 13: Results for the gauge d t a = —2 aK + ft dj a, dtft =? 
f.f * + ft dj ft — 7] ft at time t — 30 M. The results are nearly 
identical to those depicted in Fig. 10 

Numerically ft has proven to be very well behaved so f ' 1 
evidently acts as an effective “shock absorber” , at the 
expense of not vanishing entirely. 

V. CONCLUSIONS 

We have investigated gauge conditions that are appro- 
priate for the moving puncture approach to black hole 
simulations. Potential hazards from zero-speed or slow- 
speed modes have been identified, a methodology for an- 
alytically exploring various gauge choices has been pre- 
sented, and a gauge free of slow-speed modes has been 
recommended. Several gauges were also studied numer- 
ically and two gauges distinguished themselves as par- 
ticularly well-adapted to smooth black hole motion. In 
both cases we have suggested a simplification of the shift 
evolution equations, as well as explored the possibility of 
eliminating a traditional damping term to better realize 
a “Gamma-freezing” condition. We found that elimating 
this damping term on the shift does indeed minimize P 1 
and yield smooth evolutions, more so with the addition 
of our recommended “shifting-shift” terms, and warrants 
further experimentation. We intend to investigate this 
rj -* 0 limit more thoroughly in future work. 
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